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We reconstruct B° — > D~ D + decays using a sample of 449 x 10 6 BB pairs recorded by the 
Belle experiment, and measure the branching fraction to be B(B° — > D + ) = [7.5 ± 0.2(stat) ± 
0.8(syst)±0.8(£>'s)] x 10~ 3 . A search for the related decay B° — > D~ is also performed. Since we 
observe no statistically significant signal an upper limit on the branching fraction is set at 3.6 x 10~ 5 
(90% C.L.). 

PACS numbers: 13.25.Hw, 14.40.Nd 



I. INTRODUCTION 

Several decay modes of B mesons with a in the final 
state have been measured at the B-factories. The ampli- 
tudes governing these decays are interesting because none 
of the constituent flavors of the are present in the ini- 
tial state. For example, the decays B°-> DfK~ [l|, Q 
and B° -> D* (2317)+ K~ 0], observed with branching 
fractions in the range 10~ 5 — 10~ 4 , can proceed via a 
bd — ► cu VF-exchange diagram. Here we study the re- 
lated decays B° -> D+Dj and B° -> D~ D+ . The for- 
mer proceeds via Cabibbo-suppressed t-F-exchange and 
has not yet been observed; theoretical calculationspre- 
dict a branching fraction ranging from ~8 x 10~ 5 [J] up 
to ~3 x 10~ 4 [5]. The latter of the two above decays 
proceeds via a Cabibbo-favored tree diagram; the ratio 
of its branching fraction to that for B° — > D + ir~ can 
be used to test the factorization hypothesis for exclusive 
non-leptonic decays of B mesons |6|. However, previous 
measurements of B(B° -> D+D~) p], i, i, [ijj have large 
uncertainties, which limit the usefulness of this method 
at present. 

In this paper we report an improved measurement 
of B° — > Dj D + decays and a search for B° — > 
D+Dj decays with the Belle detector [TO] at the KEKB 
asymmetric-energy e + e~ collider Charge conjugate 
modes are implied throughout this paper. The results 
are based on a 414 fb _1 data sample collected at the 



center-of-mass (CM) energy of the T(45) resonance, cor- 
responding to (449.3 ± 5.7) x 10 6 BB pairs. We assume 
equal production of B°B° and B + B~ pairs. To study 
backgrounds, we use a Monte Carlo (MC) simulated sam- 
ple [13[ of T(45) — * BB events and continuum events, 
e + e~ —>■ qq (q = u, d, s and c quarks). 



The Belle detector is a large-solid-angle magnetic 
spectrometer that consists of a multi-layer silicon ver- 
tex detector (SVD), a 50-layer central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters 
(ACC), a barrel- like arrangement of time-of- flight scintil- 
lation counters (TOF), and an electromagnetic calorime- 
ter (ECL) comprised of CsI(Tl) crystals located inside a 
superconducting solenoid coil that provides a 1.5 T mag- 
netic field. An iron flux-return located outside of the 
coil is instrumented to detect K\ mesons and to iden- 
tify muons (KLM). The detector is described in detail 
in Ref. [111 ]. Two different inner detector configurations 
were used. For the first 152 million BB pairs, a 2.0 cm ra- 
dius beampipe and a 3-layer silicon vertex detector were 
used; for the latter 297 million BB pairs, a 1.5 cm ra- 
dius beampipe, a 4-layer silicon detector and a small-cell 
inner drift chamber were used 1141. 
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II. RECONSTRUCTION 

Charged tracks are selected with loose requirements on 
their impact parameters relative to the interaction point 
(IP) and the transverse momentum of the tracks. For 
charged particle identification (PID) we combine infor- 
mation from the CDC, TOF and ACC counters into a 
likelihood ratio C(K)/(C{K) + £(tt)) 0. A selection 
imposed on this ratio results in a typical kaon (pion) 
identification efficiency ranging from 92% to 97% (94% 
to 98%) for various decay modes, while 2% to 15% (4% 
to 8%) of kaon (pion) candidates are misidentified pions 
(kaons). 

We use the D~ -> 07r~ , K*°K~ and K^R- modes to 
reconstruct Dj mesons and D + — > K + K~ir + , K~tt + it + , 
and Kgir + for the D + mesons, where the <fi, K*° and 
Kg decay to K + K~, K + tt~ and 7r~7r + , respectively. 
Combinations of oppositely-charged kaons with (m^ — 
M K + K - | < 20 McV/c 2 and of oppositely-charged kaons 
and pions with \m K *o — M K +^-\ < 85 MeV/c 2 , orig- 
inating from a common vertex, are retained as 4> and 
K* a candidates, where and m K ,a are the nominal 
masses of the two mesons [lq |. Neutral kaons (K%) are re ~ 
constructed using pairs of oppositely-charged tracks that 
have an invariant mass within 30 MeV/c 2 of the nomi- 
nal K° mass, and originate from a common vertex, dis- 
placed from the IP. All Dr s ) candidates with invariant 
masses within a 4 a (4.5 a) interval around the nominal 
D s (D) mass arc considered for further analysis, where 
D s (D) signal resolutions (er) range from 3.6 MeV/c 2 to 
4.2 MeV/c 2 (3.7 MeV/c 2 to 4.1 MeV/c 2 ). A decay ver- 
tex fit with a mass constraint is applied to the selected 
D( s ^ candidates to improve their momentum resolution. 
For the decay B° ^ D+Dj we also add an additional 
constraint on the value of the cosine of a helicity angle, 
| cos fc | > 0.05 (0.25) for the D~ -> fa- (K* q R-) de- 
cay mode, where Oh is defined as the angle between the 
direction of the D~ and the K + originating from the 
vector-meson (<fi or K*°) in the vector-meson rest frame. 
The distribution in cos Oh is expected to be proportional 
to cos 2 9h for the signal and uniform for the combinatorial 
background. 

Pairs of D~ and Df , meson candidates are com- 
s \ s ) 

bined to form B° meson candidates. These are iden- 
tified by their CM energy difference, AE = E^ M — 
£™ m , and the beam-energy constrained mass, M bc = 

VWcaJ 2 " (pg M ) 2 , where E™ m = is the CM 

beam energy and E^ M and p§ M are the reconstructed 
energy and momentum of the B meson candidate in the 
CM frame. The signaPregion is 5.272 GeV/c 2 < M hc < 
5.285 GeV/c 2 for the B° -> D~ D+ , and 5.274 GeV/c 2 
<^ M bc < 5.284 GeV/c 2 and \AE\ < 0.013 GeV for the 
B° -> DfD- decays. 

To suppress the large combinatorial background domi- 
nated by the two-jet-like e + e~ — > qq continuum process, 
variables characterizing the event topology are used. We 
require the ratio of the second to zeroth Fox- Wolfram mo- 
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FIG. 1: AE distribution for reconstructed ~B° -> D~ D + 
events in the M^ c signal region. The curve shows the re- 
sult of the fit. The normalized distribution for the events in 
the sidebands of both D 3 and D invariant masses is shown as 
the hatched histogram. 

ments [17[, i?2 < 0.3 and the thrust value of the event, 
T < 0.8. Simulation shows that this selection retains 
more than 95% of BB events and rejects about 55% of 
cc events and 65% of uu, dd and ss events. 

The above selection criteria and signal regions are 
determined by maximizing the figure of merit (FoM), 
S/ y/S + B, where S and B are the numbers of signal and 
background events determined from MC. For optimiza- 
tion of the FoM we assume B(B° -> D+D~) = 2 xl0~ 4 . 

The fraction of events with more than one B — > 
D-D+ (B° -> D+D~) candidate is 4.9% (2.8%). As 
the best candidate we select the one with the mini- 
mal x 2 — X 2 (-D7) + X 2 (-C(s)) value, where x 2 (Z?j) and 
% 2 (D^) are % 2 's of the mass-constrained vertex fit. 

III. b°^d;d+ decays 

The AE distribution of events in the Mb c signal region, 
obtained after applying all selection criteria described 
above is shown in Fig. [T] Apart from the signal peak 
at AE = 0, contributions from two other specific de- 
cay modes were identified using the MC: B° — > D*~D + 
and B° — * D~ D* + . These events cluster around AE 

0.1 fi GeV and —0.10 GeV due to the unreconstructed 
7T° or 7 from the D* s ^ meson. 

The AE distribution is described by two Gaussians 
with the same mean for the signal, two Gaussians for the 
B° — > D*~ D + , D~ D* + background events, and a linear 
function for the rest of the background. The normal- 
izations, positions and widths of the Gaussians are free 
parameters of the binned likelihood fit. The solid line 
in Fig. [1] shows the result of the fit. The positions and 
widths of the B° — > D*~ D + , D~ D* + background com- 
ponents agree with the values expected from the MC. 
In addition, we perform separate fits to the AE distri- 
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butions for each D s decay mode using the same function 
with the widths and means of all four Gaussian functions 
fixed to the values obtained by the overall AE fit. 

We use events in the D s and D meson invariant mass 
sidebands in order to check for peaking backgrounds. For 
this check the masses of D s and D candidates are not 
constrained to their nominal masses. The D(D a ) invari- 
ant mass sidebands are ±200 MeV/c 2 intervals around 
D{D S ) nominal mass, excluding the D{D S ) signal region. 
Due to common final states used to reconstruct D and 
D s candidates we exclude the D S (D) signal regions and 
a ±27 MeV/c 2 D* + mass region from D(D S ) sidebands. 
The AE and M^ c distributions obtained by simultane- 
ously using events in the sidebands of both the D and D s 
mesons are in agreement with the observed combinatorial 
background under the B — > Dj D + signal. A significant 
signal is present only in the D s sideband, for D s 's recon- 
structed in the D~ — > K*°K~ decay mode. This is due 
to the three-body B° — > D + K*°K~ decay, reported in 
Ref . . The fraction of these events in the signal peak 
was evaluated by fitting the AE distribution in the D s 
sideband. We observe no peaking background when using 
the D mass sideband. The signal in Fig. [T] also includes 
contributions from Df — > K + K~tt + , / O (980)7r + and 
Kq(143Q)°K + , which all have a common final state, as 
well as a small contribution (0.4%) from D + — ► K~it + it + 
decays, where one of the ir + decays in-flight to a /i + and 
i/fj, and the /i + is misidentified as the n + . We evaluate 
these fractions using simulated events. The contribution 
of these decays is around five times larger than the con- 
tribution of B° — > D + K*°K~ decays. We take into ac- 
count the relative contributions of individual D s and D 
decay modes and determine the overall fraction of peak- 
ing background events (r) to be (11.3±2.6)%. The uncer- 
tainty includes the statistical uncertainty in D s sideband 
fits, non-uniformity of M(K*°K~) in B° -> D+K*°K~ 
decays, limited MC statistics and uncertainties in the 
corresponding branching fractions (16J. 

The signal yield for B° -> DjD + is thus N = (1 - 
rJiVpeajj = 2230 ± 56(stat), where N pea k is the number of 
events in the signal peak obtained from the fit to the AE 
distribution (Fig.[T|). 



IV. B°-» DECAYS 

The AE distribution for B° -> D+D- decays ob- 
tained after applying all selection criteria described above 
is shown in Fig.^a). The expected width of the narrower 
signal Gaussian, which describes 82% of the events, is 
5.8 MeV. This value is obtained from the MC sample 
and rescaled by a factor obtained after a comparison of 
parameters from B° — ► D~ D + data and MC samples. 
The AE signal region includes around 89% of the signal. 

While the AE distribution of the combinatorial back- 
ground is well described by a first order polynomial, 
there is a significant cross- feed contribution from B° — > 
D~D+, D*~D + , and D~ D*+ decays, where the D+ dc- 
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FIG. 2: (a) AE distribution for the 73° -> D+D~ decay 
mode. Two vertical dashed lines show the interval excluded 
from the fit, as described in the text, and two dotted lines 
show the AE signal region, (b) AE distribution for recon- 
structed events obtained by inverting the kaon identification 
requirements in data and in the MC sample. 



cays into a K~n + ir + or K°ir + final state and one of the 
pions is misidentified as a kaon. Figure [2fb) shows the 
AE distribution of these cross-feed events, as obtained 
in both data and MC samples by selecting one of the 
kaon tracks in the D s decay chain with a pion PID re- 
quirement. Events peaking around 0.075 GeV are due to 
B° — > D~D + decays, while the events clustering around 
-0.1 GeV are due to B° -> D*~D + and B° -> Dj D*+ 
decays without a reconstructed ir° or a photon. The AE 
distribution of cross-feed events is described by the sum 
of two Gaussian functions and a constant. The solid line 
in Fig. Hfb) shows the result of the fit. The widths and 
means of the two Gaussian functions are statistically con- 
sistent with the values obtained from MC. 

The expected number of background events populat- 
ing the AE signal region is determined by a binned like- 
lihood fit to the AE distribution sidebands (|A_E| > 
24 MeV region indicated by the two vertical dashed 
lines in Fig. [5Ja)). While normalizations are free pa- 
rameters of the fit, the widths and means of the two 
Gaussian functions are fixed to the values obtained from 
a fit to the AE distribution of the misidentified data 
(Fig. IDJb)). The fit result is then integrated across the 
AE signal region (indicated by the two dotted lines in 
Fig. EKa)) to obtain the number of background events, 
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TABLE I: Sources of systematic uncertainty in B(B° 
DJ D + ) and B(B° — > Df D~) measurements. 



Systematics 




"[%] fl°->£>+£>7[%] 


£Ts of D a and D mesons 


10.1 


18.6 


Tracking 


6.0 


6.0 


PIDiK^/^/Kl t 


7.1 


8.3 


MC statistics 


1.7 


3.9 


Signal window 


1.0 


2.1 


Signal fraction (1 — r) 


2.9 




Fitting model 


1.9 


included in b 


N(BB) 


1.3 


1.3 


Total 


14.5 


21.7 



b = 6.7±0.8(stat)±0.5(syst), where the systematic error 
is evaluated by varying values of the fixed fit parame- 
ters by one standard deviation. Since only three events 
are observed in the AE signal region, the result for b 
indicates that there is no statistically significant signal 
present in this AE interval. Thus the expected 3.5% tail 
of the signal, which might populate the fitted region (pa- 
rameterized as background only), can be safely neglected. 

The average efficiency of the selection criteria 
e(D s D s ) = J2 h3 e^B{D sl )B(D SJ ) = (2.45 ±0.46) x 10~ 4 
is evaluated from MC, where the intermediate branch- 
ing fractions B(P S — > 07r) and B(D S — > K*°K) are taken 
from Ref. 0, and B(D S -> K° S K) is taken from Ref. [H|. 

To check for a possible peaking background we use 
events in the D s mass sidebands. No peaking structures 
are observed in any of the M^ C -AE distributions. 



V. RESULTS 

We consider several sources of systematic uncertainty 
as listed in Table |U The largest contribution arises from 
an imprecise knowledge of the intermediate branching 
fractions of D s and D mesons and amounts to ±10.1% 
(±18.6%) for B° -> D-D+ (5° -> D+Dj) decay 
mode 16, [H. A 1% relative error for each of the charged 
tracks used in the reconstructed final states is assigned 
due to the uncertainty in tracking efficiency determined 
using partially reconstructed D*'s. The particle identifi- 
cation efficiency has a relative uncertainty of 1.4% per 
charged kaon and 0.8% per charged pion, determined 
from D*+ £>°7r+, P° -> K~-k + decays. The rela- 
tive error for each reconstructed Kg in the final state is 
4.5%. A 1.7% (3.9%) uncertainty is due to the limited 
statistics of the MC sample used for the efficiency cal- 
culation. Since the efficiency is evaluated for the signal 



region, we assign an additional 1% (2.1%) uncertainty 
due to the small possible difference in the signal resolu- 
tion between data and MC samples. A 2.9% uncertainty 
is due to the imprecise knowledge of the fraction of true 
signal events, (1 — r), in the data sample. Systematic 
uncertainty arising from the description of the AE dis- 
tribution is evaluated by comparing the known number 
of reconstructed B° -> D~ D + events in the simulated 
sample with the fitted yield and is found to be 1.9%. Fi- 
nally, the uncertainty in the number of BB events (1.3%) 
is taken into account. The sum in quadrature of the in- 
dividual contributions gives a systematic error of 14.5% 
for a B(B° -> D~D + ) and 21.7% for a B(B° -> D+Dj) 
measurement, excluding the uncertainty due to the fit- 
ting model included in b. 

The number of signal B° — > PjP + events, Nd b d, is 
converted into a branching fraction using the MC effi- 
ciency e(D s D) and the number of BB events. The mea- 
sured branching fraction is given in Tabic HH 

We use the world average of B(B° — * D + ir~) and 
calculate the ratio 



pex. 
H D 3 /ir 



B(B° -> D-P^ 
B(B l 



>P + TT- 



= 2.65 ±0.42. 



(1) 



Before comparing this result to the numerical predic- 
tion of Rd s /tt given in Ref. Q — in which the calcula- 
tion is performed in the generalized factorization scheme 
and includes penguin effects — we rescale it by a fac- 
tor {fo^/fof) 2 , where f™ w is the average value of D s 
meson decay constant given in Refs. [TH, [2(| and /£, ld is 
the value used in the original calculation. The expected 
value is R^' j n — 3.12 ±0.35, where the uncertainty orig- 
inates from the dependence on the decay constant fo s 
and form-factors, the former being the main source. The 



ratio R°n i 



k/ R d' b /k = °- 85 ± 0.13(ex.) ± 0.09(th.) is con- 
sistent with unity. If one does not include the penguin 
contributions [|[ to the amplitude for B° — > P~P + de- 
cay, the above ratio would be 0.61±0.10(ex.)±0.07(th.). 

We observe no statistically significant signal in the 
B° — > PfPj decay mode. The central value for 
the measured branching fraction is [—3.4 ± 1.6(stat) ± 
0.6(syst) ± 0.6(B's)] x 10~ 5 . We infer an upper limit 
on the B(B° — > PfP~) from the total measured num- 
ber of reconstructed events and the number of back- 
ground events in the AE signal region (no = 3 and 
b = 6.7±0.9, respectively), and the measured sensitivity, 
So = N BO g B -e(P s P s ) = (110±24)xl0 3 . The latter error 
includes all systematic uncertainties given in Table [H To 
estimate the upper limit we use Bayes's theorem with a 
flat-prior for the si gna l following the prescription in (Sec- 
tion 32.3.1 in Ref. [Hi): 



J 



p(B\n ,b,a b , S ,a s ) = 



J^o X!^o L(n a \B, fJ.b)G(S\S , (j s )tt{B, n b \b, a b )d/j, b dS 
SZ, IZ, I-oc L(noW> tiW'\S , as)n(B', y! b \b, a b )dB'd^' b dS' ' 



(2) 
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TABLE II: Results on the fitted numbers of events in the signal peak and branching fractions for B —> DJ D decay mode. 
The peaking background fraction, r, is given for each D a decay mode in the second column. The efficiencies include intermediate 
branching fractions (e(D s D) = £jB{D a )B(Dj)), which are taken for all three D and D a — > (pit and K*°K decay modes from 
Ref. [ID], and that for D s — > K%K is taken from Ref. In the last column the dominant sources of systematic uncertainty, 
the D s branching fractions, Bd s x, are factored out. Quoted uncertainties for B are statistical, systematic and uncertainty due 
to the imprecise knowledge of intermediate branching fractions, respectively. 



Mode 


N peak 


r[%] 


e [10" 4 ] 


B [10" 3 ] 


B ■ B DsX [10" 4 ] 


D- -> 071- , (j>^ K+K~ 


1112 ± 35 


12.9 ± 4.5 


2.8 ±0.4 


7.8 ± 0.2 ± 0.9 ± 1.0 


1.68 ±0.05 ±0.19 ±0.06 


D~ -> K*°K~, K*° -» K+iv- 


961 ± 33 


14.5 ± 4.3 


2.5 ±0.5 


7.3 ±0.3 ±0.8 ±1.5 


1.83 ±0.06 ±0.21 ±0.07 




441 ± 22 


0.4 ± 2.2 


1.3 ±0.1 


7.3 ±0.4 ±0.9 ±0.6 


0.76 ±0.04 ±0.09 ±0.03 


Combined 


2514 ± 64 


11.3 ± 2.6 


6.6 ±0.7 


7.5 ±0.2 ±0.8 ±0.8 





The number of observed events no is Poisson distributed 
around the sum of ^ s and ^b- L(rio\fi s , = 1/no! (Hs + 
Mb)"° e ~U*'+M>) ) where n s and in, are the expected number 
of signal and background events, respectively. In partic- 
ular ij, s can be written as fi s = B ■ S, where B and S 
are true values of B(B° — > D+D~) and the sensitiv- 
ity N B ago ■ e{D s D s ), respectively. The true value of S 
can only take non-negative values and is Gaussian dis- 
tributed around Sq with variance erg. Hence Q(S\So, as) 
is a Gaussian function with a cut-off for S < 0. The prior 
probability density n(B, /xj,|6, <7b) is assumed to be factor- 
izable, n(B, Hb\b, <Jb) = P(B)Q(nb\b, cr&). For P(B) we use 
a flat-prior, and G((J>b\b,ab) is again a Gaussian function 
centered at 6, with a width of Ob and with a cut-off for 
b < 0. 

Integrating out the nuisance parameters S and /ib we 
obtain the posterior p(B\hq, b, <jb, So, erg), which already 
takes into account the statistical error on 6, the system- 
atic error due to the parameterization of AE distribution 
in the fit, and systematic uncertainties on the efficiency 
and on the number of BB pairs. The 90% C. L. upper 
limit on B(B° — * D+D~) following from this posterior 
is found to be 

B(B° -» D^DJ) < 3.6 x 1CT 5 at 90% C.L.. 

VI. CONCLUSIONS 

In conclusion, we have measured the branching fraction 
for B — > DjD + decays. The measured value is B(B — > 



DjD+) = [7.5±0.2(stat)±0.8(syst)±0.8(#'s)] x 10~ 3 , 
which represents a large improvement in accuracy as 
compared to previous measurements 0, H, H, E3] ■ Com- 
bining this result with the world average for B(B° — ► 
D-JT+) we obtain the ratio RdI^/^d,/^ = 0-85 ± 
0.13(ex.)±0.09(th.). With present experimental and the- 
oretical uncertainties, the results are consistent with the 
factorization hypothesis for non-leptonic exclusive decays 
of B mesons. If one does not include the penguin con- 
tributions @ to the amplitude for B° — > D~ D + decay, 
the above ratio is not consistent with unity. For B° — > 
Df D~ decays we found no statistically significant signal. 
We set an upper limit of B(B° -> D+Dj) < 3.6 x 10~ 5 
at 90% C.L. This result puts even more stringent limits 
on 8(B° — > D+D~) than the recent measurement by 
the BaBar collaboration [2l[, severely challenges recent 
theoretical estimates in Refs. 0, H| and implies that the 
weak annihilation contributions in decay modes with two 
charmed mesons are small, as suggested in Ref. [221 ]. 
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